Reentrant behavior of superconducting alloys 
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A dirty BCS superconductor with magnetic impurities is studied. Asymptotic solution of the ther- 
modynamics of such superconductor with spin 1/2 and 7/2 magnetic impurities, is found. To this 
end, the system's free energy f{H,P) is bounded from above and below by mean-field type bounds, 
which are shown to coalesce almost exactly in the thermodynamic limit, provided the impurity con- 
centration is sufficiently small. The resulting mean-field equations for the gap A and a parameter 
u, characterizing the impurity subsystem, are solved and the solution minimizing / is found for var- 
ious values of magnetic coupling constant g and impurity concentration x. The phase diagrams of 
the system are depicted with five distinct phases: the normal phase, unperturbed superconducting 
phase, perturbed superconducting phase with nonzero gap in the excitation spectrum, perturbed 
gapless superconducting phase and impurity phase with completely suppressed superconductivity. 
Furthermore, evidence of reentrant superconductivity and Jaccarino-Peter compensation is found. 
The credibility of the theory is verified by testing the dependence of the superconducting transition 
temperature Tc on a;. Very good quantitative agreement with experimental data is obtained for sev- 
eral alloys: (Lai_iCei:)Al2, (Lai_a;Gda;)Al2 and (Lao.8-a:Yo.2o)Cei. The theory presented improves 
earlier developments in this field. 

PACS numbers: 74.20.-z, 74.25. Bt, 74.25.Dw, 74.70. Ad 



I. INTRODUCTION 

It is well known that doping may substantially change 
the properties of a superconductor. The superconducting 
transition temperature T^ in most dirty superconductors 
decreases with impurity concentration c^. However, al- 
loys of zirconium with iron, cobalt or nickel have higher 
transition temperatures than pure zirconium^. Similar 
behavior was observed in titanium doped by chromium, 
manganese, iron and cobalt^. One of the properties, 
which do not change under doping, is the order of the 
transition in the absence of external magnetic field^i^. 

Early attempts to explain lowered T^ in the presence 
of magnetic impurities were founded on perturbation 
theory. NakamuraS. and Suhl et alJ- explained this ef- 
fect by treating the s-d interaction V^-tj^ as an additive 
term in the total Hamiltonian, which perturbs a BCS 
superconductor^. However, their theory predicts a first 
order phase transition to the superconducting state in 
zero magnetic field. 

The well known Abrikosov-Gor'kov theorjiiS, (AG) of 
dirty superconductors explains the strong decrease in 
Tc due to magnetic impurities and also predicts "gap- 
less superconductivity" , confirmed experimentally by 
Reif and Woolf^^. Disagreement with this approach is 
observed in a number of Kondo superconductors, e.g. 
Lal_:,Ce^Al2l^ LaCe and LaGd^^ and PbCe and InCei^. 

Such compunds manifest reentrant superconductivity 
(RSC) which is due to competition between the Kondo 
effect and superconductivity. The RSC effect was pre- 
dicted theoretically by Miiller-Hartmann and Zittartz— 
(MHZ) and is expected to occur when the superconduct- 
ing transition temperature of the host compound Tco is 
much larger than the Kondo temperature Tk. In such 
the case the equation for Tc(c) has three solutions: Td, 



Tc2 and Tcs for certain values of impurity concentra- 
tion. As the temperature T is lowered the alloy becomes 
superconducting at Td. When the value of T is com- 
parable to Tk the pair-breaking effect of impurities in- 
creases and superconductivity is suppressed at Tc2- The 
alloy re-enters the superconducting state at Tc3, where 
the pair-breaking passes through its maximum value. 
Such reentrant behavior was observed experimentally in 
La0.7915Ce0.0085Y0.20 with Td = 0.55 K, Tc2 = 0.27 K 
andTc3 = 0.05Ki6. 

According to Miiller-Hartmann and Zittartz theory, 
superconductivity is never completely suppressed for any 
value of impurity concentration if Tk/Tco ^ 1. This 
statement disagrees with experiment, which shows that 
disappearance of superconductivity above a critical value 
of c is possible. Furthermore, significant deviations from 
this approach were observed for (La, Ce)In3_3;Sna;^^. 

First experimental observations of reentrant 
superconductivit y^^'^^ revealed no evidence of Tc3, 
corresponding to second phase transition, reintroducing 
superconductivity. For this reason, MHZ theory was 
reformulated (e.g. Refs. l20l - |23l ) and also other proposals 
for a theory of Tc(c) were given. In particular Jarrell, 
performed Monte Carlo simulations of the superconduct- 
ing transition temperature in terms of Eliashberg-Migdal 
perturbation theory—. These calculations raised doubts 
about the existence of Tc3, contrary to experiments 
accomplished by Winzer— , confirming the presence of a 
transition back to the superconducting state. 

The effect of magnetic impurities on superconductivity 
is still under debate. Recently, Barzykin and Gor'kov^— 
studied s-wave superconductivity in the Anderson lat- 
tice and demonstrated excellent agreement of the result- 
ing Tc{x) graphs for Cei_a;Laa;Ru3Si2 with experiment. 
Reentrant behavior of Tc may occur for untypical values 
of parameters. Kozorezov et ali^ have shown, in terms 



of the MHZ model, that trace concentrations of mag- 
netic impurities may also result in significant changes in 
nonequilibrium properties of superconductors. A com- 
prehensive review of recent developments in this field can 
be found in Ref. l27l 

Experimental studies of superconductors containing 
magnetic impurities carried out by Matthias^ revealed 
another extraordinary property, viz., the coexistence of 
magnetism and superconductivity. Till then these two 
phenomena were believed to be mutually exclusive, since 
the internal magnetic fields generated in magnetically 
ordered systems are much larger than the typical crit- 
ical fields of superconductors. The coexistence hypoth- 
esis was confirmed shortly after in the following super- 
conducting alloys: Cei-xGdiRu2^^ and Yi_2;GdxOs2^, 
although it was not found to occur in the same vol- 
ume element. The phase diagrams obtained by Wilhelm 
and Hillenbrand for Cei_2:Tba;Ru2^° also contains the 
coexistence phase. Specific-heat measurements showed 
short-range ordering in this alloy, typical of spin-glass 
systems^. 

The coexistence of superconductivity and long-range 
antiferromagnetic ordering of the rare earth R magnetic 
moments was discovered in RMoeSeg (R = Gd, Tb and 
Er)32, RRh4B4 (R = Nd, Sm and Tm)k and in RMogSg 
(R — Gd, Tb, Dy and Er)'^^. A similar overlap between 
superconductivity and ferromagnctism was observed in 
ErRh4B4.^ and HoMogSs^^. 

The phase diagrams of superconducting alloys, con- 
taining the coexistence phase, have been computed by 
several theoretists. Gor'kov and Rusinov extended the 
AG theory to include magnetic ordering. Correspon- 
dence with the phase diagrams observed experimentally 
is expected to occur for very strong spin-orbit scattering. 
Balseiro et al.— studied a BCS superconductor perturbed 
by magnetic impurities interacting via a nearest neigh- 
bour Heisenberg potential. The resulting phase diagrams 
comply qualitatively with experiment. 

Theories of dirty superconductors contribute signif- 
icantly to our understanding of the superconductivity 
phenomenon. One can expect that further investiga- 
tions of superconducting alloys will explain the micro- 
scopic mechanism of unconventional superconductivity 
displayed by some materials, e.g. high temperature 
superconductor^ and heavy fermion compoundsSS. 

This issue, as well as some shortcomings of the models 
presented above, motivate the present work. We investi- 
gate a BCS Hamiltonian^ 



Hbcs — Hq + Vbcs , 
supplemented by a reduced s-d interaction 

M 



where 



a=l 
Ho = y ^Ck"-kcr, 



(1) 



(2) 



with fk = Ek - A*, "ko- = 
kinetic energy operator and 



'^kcr'^kff' i^ ths free fermion 



Vbcs 



-|A| 



kk' 



Gkk'aj^+a^k-"-k'-«k' 



(3) 



is the Cooper pairing potential. |A| denotes the system's 
volume and Gkk' is real, symmetric, invariant under k — > 
— k or k' — >■ — k' and nonvanishing only in a thin band 
close to the Fermi surface, viz., 

Gkk'=Gox(k)x(k'), Go>0, 
where x(k) denotes the characteristic function of the set 

?i2k2 



V = {'k: n-S <e^< H + S}, 



e^ = 



2m 



In equation ([2]), Sza denotes the spin operator of the 
magnetic ion, whereas 

a-z = ^(r^k+ - "k-) , 

k 

describes the spin operator of each conducting fermion. 
M is the number of magnetic impurities, N the number 
of host atoms. 

We assume the perturbation implemented by the lo- 
calized distinguishable magnetic impurities to be a re- 
duced long-range s-d interaction, which involves only the 
z-components of the impurity and fermion spin oper- 
ators (Eq. ^). The reason for this simplification is 
that the thermodynamics of the resulting Hamiltonian 
H — Hq + Vbcs + V admits a mean-field solution which 
improves with decreasing impurity density. Furthermore, 
this solution is thermodynamically equivalent to the one 
obtained for H with a Heisenberg type reduced s-d inter- 
action 



V„ 



N ^^ 

ko: 



[(4-ak 



-»3^ — 



^k+"k4 



_5„ 



)S. 



(4) 



■'a± 



^r\rr i lO, 



Qyi 



replacing V (Ref. [2J, Sec. 6.2.5). (Similarly, the thermo- 
dynamics of a classical superconductor can be explained 
in terms of a reduced BCS interaction, whereas a gauge- 
invariant theory of the Meissner effect requires a more 
general pairing potential.) The reduced form of Vh, ob- 
tained by rejecting the sum ^ Vkk' in the s-d interac- 

k^^k' 

tion Vs-dj is an approximation which resorts to the fact 
that spin-exchange processes, and not momentum ex- 
change processes, are primarily responsible for the Kondo 
effect caused by Ki-d^°. 

In Sections |TT] and IIIII the upper and lower bound to 
the system's free energy f{H, 0) are derived by exploiting 
the method developed in Refs. l4l| - |44l These bounds are 



shown to be almost equal if the impurity density Af/|A| 
is sufficiently small (Sec. lIVp . The mean- field equations 
for the gap A and parameter v, characterizing the im- 
purity subsystem are derived in Sees. IVHVll and solved 
in Sec. I VIII for various values of g, M, Go and 6. In 
section IVIIII the values of these parameters are adjusted 
to fit the experimental Tc{x) curves for ha,i-xCexAl2, 
Lai_a;Gd2;Al2 and Lao.8-3;Yo.2oCe2;. Quantitative agree- 
ment is found for each alloy. Furthermore, the phase 
diagrams derived in Sec. |VlT] qualitatively reproduce the 
experimental observations of the coexistence phase, reen- 
trant behavior, gapless superconductivity and Jaccarino- 
Peter compensation. 



II. LOWER BOUND TO THE FREE ENERGY 

A version of the Tindemans and Capel method^i^, 
introduced recently in order to study thermodynamics 
of the Fermi gas interacting with randomly distributed 
magnetic impurities^, will be applied in this section to 
derive a lower bound to the system's free energy. 

The full Hamiltonian of the system is 

ij(*^) = To + 1/bcs + V. (5) 

Exploiting the identity 

2 a 12/ rt\2 ,122,12 rt2 

g CTzOza = -2-9 i'^2 ~ ^za) + 29 ^ z + 2^ ^ zai 

the partition function can be written in the following 
form 



Z{M) = Trexp(-/3i/(^^)) = Trexp 



-/S^BCS 



.2p 



M 



a=l 

The Gaussian integral 

1 ' 



exp la 



/2^ 



Sza) — Oy, 



exp ( —\v^ 



S: 



(6) 



\j2av ) Av 



where -ffecs = ^o + Vbcs ■ In order to linearize quadratic 
terms in u^. hi the exponent of the integrand on the r.h.s. 
one exploits the inequality 

iiV-igVf = \N-^ [go, - N^f + g^a, - \Nif 
> gijcT, -\NTf, 

where 77 is an arbitrary function of z^i , . . . , vm ■ This yields 
an upper bound to the partition function 

Z{M)<[^] Trexp [-/3Hbcs] 



\2n J 



M 



' a = l 



Y\^ di/Q exp -^N ^g'^jSSl^ - gPuaS^a 



\Np{yl-if)+f3g{v^~-j^)a, 



(9) 



Following Pearce and Thompson^, let us now add and 
subtract the term ^Nfix^^v^^^ with a; > 1 in the expo- 
nent of the integrand. The resulting expression is next 
split into two factors and one of them is replaced by its 
maximum with respect to {i'q}: 

M/2 



27r ) {u^} 

M 



Z{M) < ( ^ ) maxexp(G'(z^i,...,i^Af,»7)) 



Y[di^a' exp [-iN^i^Kl - x'^)] 

-oc ^, 

maxexp(G(i/i,...,i/M,'7)) {^ - x~^) 



-M/2 



(10) 



where 



G{vi,...,VM,ri) = In Trexp 



-/3iJi 



BCS 



M 



(11) 



and the commutation relations 



Vbcs, ^n^,"- 



Vbcs, ^n\^" 



0, (7) 



satisfied for Gkk' of the form Gkk' = Gox(k)x(k'), allow 
to separate the electron and impurity spin operators in 
®, 

^(M)=(^) Tr|exp[-/3i?BCs] 



M 

1 j Ava exp 



-^fiVa + PgVa {o-z ~ Sza) 



\PN-^g^al - \pN-^g^Sl^ 



(8) 



\N-^g''Sl^ - giy^Sza + hNif 



Inequality (|T0|) yields the relevant lower bound to the free 

energy|A|/(iJ,/3) = -/3-ilnZ, 

|A|/(iJ,/3) > - min/?-iG(z.i, . . .,z/M,77) 

{u^} (12) 

+ \l3-Hl\n{l-x-^). 

The function 77 will be now chosen as the solution of the 
equation: 



dG_ 

drj 



= 0, 



which according to pT|) reduces to 

j^^gN-^{az)r. 



(13) 



(14) 



where 



IM 



Hbcs ~ g'^ii^a -ri)crz, 



(15) 



/,(i!f) = g.y^5.. + iiv-vE^- 



imp 






(25) 



and 



{A)h 



Tr(Aexp[-/3i7]) 
Trexp[-/3H] ' 



(16) 



The sohition 'q{vi,--- ,vm) of equation (|14|1 is unique. 
The proof can be found in Refs. l42H43l . For this choice of 
77, the condition for the minimum in (J12p 



dG 



0, A = 1,...,M, 



simphfies to 



_i _ g Tr5';,A exp [- a^/^g^g*^;, - PgS^xvx] 
"^ "'^ N TTexp[-^(3g^Sj^-pgS,xi^x] (17) 
+ 5^-1 (a,)^^^, A = l,...,Af. 

The form of these equations is independent of A, thus 

J^i = • • ■ = i^M = i^- (18) 

As a consequence, Eqs. (|17l) simphfy to the form 

^_i^ ^ _ g Tr^, exp [-^^g^S'', - pgS^u] 

N Trexp[-2^/3g252_^g5^^] (^g) 



where 



/i(i/, 77) = Hbcs - gM{u - ri)(T^. 



(20) 



From Eqs. pT|) and p^ one obtains a suitable lower 
bound to the system's free energy: 

\A\f{H,l3) > -mm(3-^G{iy)+^l3-Hnn{l - x-'^) ,x > 1, 

(21) 
where G(i^) denotes the function defined in (fTTj) for i/i = 
■ ■ ■ — I'M = V and 77(1^) satisfies Eq. ([Til) . 



III. UPPER BOUND TO THE FREE ENERGY 
BY BOGOLYUBOV'S METHOD 

An upper bound on /(i?, /3) can be expressed in terms 
of the hamiltonian hP'^^(y^r\)^ which is related to ij(*^) 
by the following formulae 



i7W=/iW(,.,77) + ff(^) 



R ' 



(22) 



where 



h^^'\v,n) = h + ]i^ + iMiV(i.2 - ^\ (23) 



/l = i?BCS + KCTz, 



-gA/(i. - 77) (24) 



H 



(M) _ 



R 



vN 



iiV-iE[g(a,-5,„) 

OL 

Bogolyubov's inequality 

F{H^+H2)<F{H^) + {H2)h,, 
with Hi = h^^'\v,ri), yields 

ni?(''\/3)<n/i(*^^(^,^/),/3) 



(26) 



(27) 



(28) 



The inequality Ty{pA^) < (Tr(pA))^ , vahd for any 
bounded self-adjoint operator A and density matrix p, 
allows to bound from above the second term in (l28l) 



From Eq. ([T4|) one obtains 

{{ga, - VN)'),^,,,, < v'N' ~ ^v'N^ + v'N' = 0, 
which yields the relevant upper bound to the free energy 

F(7j(*^),/3) < ^(/7(*^)(7.,r;),/3) = -r'G{,,). (29) 



IV. THERMODYNAMIC EQUIVALENCE OF 

^(M) AND h^"'^ 



By passing to the limit |A| -^ oo in Eqs. (PTjl . (1^ . and 
subsequently x— >-l-)-e, e>Oin Eq. pT|) . one finds that 
the upper and lower bound on /(i/, /3) coalesce almost 
exactly, if the impurity density is sufficiently small, viz., 

lim d-lim/(i/(*^),/?)« lim d-\im f{h(^'\j,^,7j^), (3), 

|A|— >CXD |A|— >oo 

(30) 
where d-lim denotes the limit of small c and i^^, r]^ are 
the minimizing solutions of the equations 



iy^gN-^{a,)~^ 



9 



TrS,exp[-^/3g^S^,-l3gS,iy] 



N Trcxp[-2^/3g25|-/?g5,7.] 



77 = gAf-M^.>^ 



(31) 



(32) 



On these grounds we shall assume that the thermody- 
namics of the original system, characterized by H^^^> is 
equivalent, under these restrictions, to that of h^^'^\ 



The two equations ((3T|) . ((32|) reduce to a single one for 
1/ if 5 > 0. The general form of Eqs. §1^1, ^ is 



Let g > 0, then /2 > 0. Furthermore, 

■r] = v- !2{v), 
which yields the equation for v. 

i' = /l(/2M)+/2M, 

where according to Eqs. ([55)1 . (IMl) : 

/i(^) = (MiV/3)-i^lnTrexp[-/3/i(i.,0)], 



(33) 
(34) 
(35) 
(36) 

(37) 



h{v) = (iV/3)-i|- lnTrexp[-/3/i(;^p]. (38) 



V. MEAN-FIELD DESCRIPTION OF h 



The form of the Hamiltonian h p4|) is analogous to 



HbcsCH) = Hq + l^BCS - MbHct^ 



(39) 



describing a system of electrons with attractive BCS in- 
teraction in the presence of an external magnetic field 
Ji {fiB denotes the Bohr magneton). The explicit 
form of the system's free energy f{h'-^^'>{h',r]), (3) can be 
therefore derived by exploiting the Bogolubov-Valatin 
transformation4^'47 g^^^^ ^^le method developed in Ref. 44 
for HBCsi'H). The Bogolubov-Valatin transformation, 



yields 



k 



a'_^_a_ 



(40) 



(41) 



One expects the energies of the states k+, — k— to be 
different, since the products q;J^.^q;^.^, aLk-Q^-k- appear 
with the opposite signs in Eq. (|4ip . Thus, the trial equi- 
librium density matrix, approximating Z^^ exp[~/3h], 
has the form (see Ref. l44r ) 



Po 



exp[-^iJo] 



Trexp[-/?iJo] 
is characterized by the following Hamiltonian 



(42) 



Ha = 5I(^kia 



k+"k^ 



Eu-?a 



k2"_k- 



-)+^o, (43) 



where Eq denotes the ground state energy. According to 
(1231) and ([50)1 the system's free energy can be decomposed 
into three summands 



where 



F = Fei+ F™p + \MN {v^ - if) , (44) 



F»mp = ~r^ lnTrexp[-/3/i,„p], (45) 



is the free energy of impurities and 

F,i = Uei - TSeu (46) 

denotes free energy of electrons, with 

Uei = Tr(/ipo) (47) 

and 

5e, = -(r/3)-iTr(polnpo)- (48) 

For po defined by Eq. (ji^ . 

S,i = - (/3T)-i Y. [/ki In/ki + (1 - ./ki) ln(l - /ki) 



+ /k2ln/k2 + (l-/k2)ln(l-/k2) 



(49) 



where 



fki 



exp(-/3£:ki) 
l + exp(-/?£;ki)^ 



rk2 



exp(-/3£'k2) 



1 + eM-pE^2) 



(50) 



The projectors on the ground and excited BCS states 



Pok = \BCS)^ k {BCS\ , 
Pk+ = 4+|S^^>kk(5^'5|«k+, 
-k-l^^'5>kk(-5C5|a_k-, 

|5G^>kk(5C'5|"-k-«k+, 



"k+"-k- 



-fk+,-k- 

allow to rewrite Eq. p2)) in the following form 

PO =n((l - •^kl)(l - /k2)P0k + /kl(l - /k2)Pk+ 



+ /k2(l - /kl)P-k- + ./kl./k2Pk+,-k- 



(51) 



Thus, one obtains 

Uel - Y. [•^kl(ek + Z^) + /k2(a - k) + 2ek«k(l - /kl 
k 

- /k2)| - |A|"^^Gkk'MkWkUk'Wk'(l - /kl - /k2) 

(52) 



kk' 

X (1 -/k'l -/k'2), 



which yields 
Pel = Uei + r^ Yl [-^ki In /kl + (1 - /ki) ln(l - /kl) 



/k2ln/k2 + (l-/k2)ln(l-/k2) 



(53) 



The parameters Wk, ^k, ^ki and i?k2, are found by min- 
imizing the free energy. The thermodynamic perturba- 
tion method of Bogolyubov et al4^ shows that for such 
choice of these parameters, the following relation holds 
up to negligible terms 

lim /(/i,/3)= hm miii/(/io,/3), (54) 

lAt^oo |At^oo{Ap} 

where h^ = Hq + kct^ and 

^p - i^r' Y. GpkWk«k(i - /kl - /k2). (55) 



Minimization of Eq. (j53p with respect to Vp yields 
= 4ep^;p(l - /pi - /p2) - 2|Ari J2 Cpk^i 

k 



dvr- 



(np-^)(l-/ki-/k2)(l-/pi-/, 



p2; 



0. 

(56) 



Then, from the normalization condition u^ + v^ = ^ and 
the demand that the density matrix po should represent 
the Fermi-Dirac distribution of free fermions for A = 0, 
one obtains 



ul 



5('+|)- «-5('-t)' ^-\/^^^- 

(57) 
The equation 



^k + K- 2^kWk + 2|A| ^^Gkk'WkWkWk'Wk' 



a/ki 

X (1 - /k'l - /k'2) + /3"'(ln/ki - ln(l - /kl)) = 0, 
combined with Eqs. (|50|) . ((55)) . yields 

^ki = a("k - «k) + « + 2AkWkWk = -Bk + K. (58) 
Analogously, for i?k2 one obtains 

■E^k2 = Ck(wk ~ ^k) - '* + 2AkUkWk ^ Ek~ K. (59) 



Using equations (l35t . (|50l) . (|57|) we get the equation for 
the gap parameter 

Ak = i|Ari^Gkk'Ak'i?k'V3(/3,i^k',ek',/2), (60) 



where 

/3(/3,i?k',a',/2) 



sinh(/3^k' 



cosh(^£;k') + cosh(5/3Af/2(j^)) ' 



(61) 

Eq. (pO)) resembles the gap equation in BCS theory. The 
convexity properties of fsi/S, Ey^> , ^k' , /2) differ in general 
from those of /bcs — ta,nh{-^(3E-k,). However, in the limit 
of extremely weak magnetic coupling, g — > 0, /a reduces 
to /bcs, viz.. 



I sinh(/3gk') 

'9=0 cosh(^£;k') + l 

_ 2sinh(i/3^k')cosh(i/3^k') 
~ 2cosh2(i;9£;k') 



(62) 



= tanh(i;3£;k') 
= /bcs- 

Equalities ^, ^, ^, §UI and dH]) now lead to the 
following expression for the free energy 

i^ei = 5] [Ck + i A^i?^ V3(/3, il^k, ^k, /2) 

k ^ (63) 

- p-'^ ln[2cosh(/3^k) -t- 2 cosh(/3K)] 



The definition (l37l) . where 77 = 0, together with Eq. ([63|) . 
yields the explicit form of the function fi{i')'. 



fiH 



eg sm]i{g /3 M 1/) 

M cosh{gl3Mv) + cosh(/3^k) ' 



(64) 



The free energy of impurities and the function /2 de- 
pend on the value of their spin. In the present work we 
examine the critical temperature T^ of alloys containing 
Ce and Gd magnetic impurities. According to Matthias 
et alii the spin of Ce ion is 1/2 and that of Gd is 7/2. 
For spin 1/2 impurities one obtains. 



imp 



-/3-ilnTrexp[-/?/ji*2] 



M 



-;3-i^lnTrexp[-/3/i, 



imp] 



(65) 



-M;3-Mn[2cosh(/35i/)] +^cg^, 



and 



/2^^W = f tanM/Sgi'). 



(66) 



Accordingly, for spin 7/2 impurities 



.(i) 



F-Z = - Mr^ In 2 exp[-24^] coah{7/3gj^) + exp[-12^] cosh(5/3(7i/) 



exp[ 



4^] cosh{3l3giy) + cosh{Pgv) 



hc9\ 



(67) 



with 



/a' V) = ]^7e^P[-24^] sinh(7/35J.) + 5exphl2^] sinh(5/35J.) + 3exp[-4^] sinh(3/3g:.) + sinh(/?g:.), (68) 



2 

where 



i? = exp[-24^] cosh{l(3gv) + exp[-12^] cosh(5;3gi/) + exp[-4^] cosh{il3gv) + cosh(Agi 

I 



VI. THE CHEMICAL POTENTIAL 

Passing from summation in Eq. (j60p over momentum 
k to integration over the single electron energies ^, one 
obtains for a sufficiently thin conduction band S: 

^ = \g^pJ ^h{P.E,£,j!f\^))A^ 5=^1/2,7/2, 

(69) 
where p denotes the density of states in V . This equation 
together, with the one for v, 



_ eg sinh(/3gM/<^V)) 

^cosh(;3gM/2^"^(i^)) + cosh(;3£;k) 



f^^i-), 



(70) 



constitutes the set of equations for A and v. The prop- 
erties of a superconductor with magnetic impurities can 
be determined by solving this set of equations, which is 
supplemented by the following condition for the chemical 
potential /i: 



E 

kcr 



Tr(nkcrPo) 



(71) 



n denoting the average number of fermions in the sys- 
tem. According to Eqs. (gHl), dSO]), dSH) and dnH), this 
condition, takes the form: 



k 



[i-|i/3(/3,i?k,a,/i'^) 



(72) 



Eq. ((7^ is analogous to the BCS equation for p, and the 
properties of /a are similar to /bcs = tanh(/3£'k/2), e.g. 
both are odd functions in ^k- The solution of Eq. (17^ is 
therefore exactly the same as in BCS theory, viz., /i — e-p- 
Numerical calculations in subsequent sections are thus 
performed under the following assumptions: 



/i = Ef, 



df 



0, 



P = PF, 



(73) 



where pp denotes the density of states at the Fermi level, 

mpp 

We have also imposed the condition 



gM f2iiy) <E^ Sit T = OK, 



(74) 



which expresses weak coupling between conduction elec- 
trons and impurities. Eq. (j69p then takes the limiting 
form 



GqPf arcsinh 



A(0) 



= 1, asT^O, 



(75) 



whereas Eq. (|70p is satisfied by i/(0) = || in this limit. 
Thus at sufhciently low temperatures T, close to OK, 
A(r) is the solution of Eq. (1751) . This is substituted into 
Eq. (j70p . which is then solved for one-fermion energies 
^ S T' by exploiting the Newton-Raphson method. The 
resulting values of i^{^) are used to obtain A(T -I- AT) 
from ([55)1 by deploying Newton- Cotes quadrature as long 
as the result is self-consistent. The resulting value of 
A(r -I- AT) is used to compute iy{^) at temperature T + 
AT. This procedure is continued until T reaches the 
specified value. To ensure numerical stability the step 
AT should be sufficiently smaU, e.g. AT = lO^^K. 

In the opposite case, viz., gM f2{v) > E]^ for T = OK 
one has 

A(0) = K0) = ^. 

In such anomalous circumstances, superconductivity is 
completely suppressed at T = K and the system is de- 
scribed only by the impurity parameter v. The interme- 
diate case, when gMf2{v) = £'k has not been examined 
numerically. 



VII. PHASE DIAGRAMS 

The system's state is characterized, according to 
Eqs. dSni and dMl), by the solution of Eqs. ^^ and 
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d70| . which minimizes F{h^'^\f3). It wiU be denoted 

by {A,„,z^m}- 

Passing from summation in Eq. (j63]) over k to inte- 



gration over the single-fermion energies ^ and exploiting 
Eq. (H^ . one obtains the foUowing expression for the free 
energy 



F(") = mm 
{A 



in|pF|A| 



5 , 



1 a2 
2 



+ Af^c-i(./(^) 



5(/2'V)+^™p + ^o(A-0)+p^52}, s = l/2, 7/2, 



de 



(76) 



where i^j„j are given by the Eqs. (1551) and (p7)) . whereas 
i?o(A = 0) denotes the ground state energy of free 
fermions. Two last terms are the contribution to the free 
energy density from one-fermion states, lying outside S. 
Equations (l69|) and ([70]) clearly possess the solution 
A = i^ = for all values of /3 > 0. At sufficiently large 
values of /3 one finds also other solutions, viz., {A 7^ 
0, t/ = 0}, {A = 0, !/ 7^ 0}, {A ^ 0, 1' 7^ 0}. Accordingly, 
we distinguish the following phases: 

— paramagnetic phase P with {Am ~ 0, J^m — 0}, 

— unperturbed superconducting state SC with 
{Am7^0,z^m-0}, 

— ferromagnetic phase F without bound Cooper pairs 
and {Ajn = 0, z^m 7^ 0}, in which impurity spins 
tend to align opposite to those of conduction 
fermions (cf. Eqs. (EH) and (125))'). 

— intermediate phase D in which superconductivity 
coexists with ferromagnetism and {Am 7^ 0, z^m 7^ 
0}. 

We define the following temperatures corresponding to 
the respective phase transitions 

— Tc, 2nd order transition SC —> P. 

— TpF, Curie temperature of 2nd order transition F 
-^ P. 

— TgcDi 1st order transition D — > SC. 

— TpDi 1st order transition D ^ F. 

— TscF, 1st order transition SC -> F. 

The set of Eqs. (p^ . ([70]) has been solved numerically 
for different values of 5, M, d and Gop-p- The parameters 
were adjusted to fit the experimental specific-heat curves 
for (Lai_2:Cea;)Al2 and LaGd and the critical field curve 
in case of ThGd^. These values are used to compute the 
phase diagrams of these alloys on the grounds of Eq. ([75| . 

The phase diagrams of (Lai_2:Ce2^)Al2, for which the 
free energy is given by Eq. (j76p with s = 1/2 and impurity 
concentration c = x/ {i — x) are depicted in Fig. [TJ The 



r 



values of M, g, (5 and GqPf are given in Table HI These di- 
agrams show the decline of D and SC phase with increas- 
ing impurity concentration. The critical temperature Tc 
is a decreasing function of impurity concentration. It fol- 
lows, therefore, that the destructive effect of impurities 
increases with c and suppresses the superconductivity if 
5 or c reaches its critical value. This complies with ex- 
perimental data, in which superconductivity is expunged 
for Ce content larger than 0.67%i^i^i^. 

The supplementary phase NC with gapless supercon- 
ductivity and Am > is present, depending on the value 
of g, as a subregion of SC, D ot F phase. The appear- 
ance of NG subregion is due to negative term —gMf2{v) 
present in Ski ((55|) . 

For sufficiently small magnetic coupling constant g, 
the NG phase for (Lai_a;Cea;)Al2 lies between D and 
F phases. Thus, the pair-breaking mechanism, which 
increases with g, induces gapless superconductivity and 
then superconductivity is suppressed for sufficiently large 
g. The system undergoes a phase transition to a ferro- 
magnetic F state (cf. Fig. [T^). 

The resulting phase diagrams also show the reen- 
trant superconductivity, which is seen on Fig. [T^ for 
g « 0.72 -s/eV. As temperature T is lowered, the system 
undergoes the first phase transition to a superconducting 
state {P — >■ SC) when T = Tc- After further cooling su- 
perconductivity is suppressed at T = Tscf and the alloy 
displays ferromagnetic properties until T > Tpo- For 
T < Tpu superconductivity reappears in the D phase, 
where it coexists with ferromagnetism. 

It is worth noting at this point that reentrant supercon- 
ductivity due to the Kondo effect was observed for the 
first time in (Lai_2;Cea;)Al2r^'i^. These measurements 
showed destruction of superconductivity below the sec- 
ond critical temperature Tc2 < Td, but the second phase 
transition, reintroducing superconductivity was not con- 
firmed. This scenario with Td = Tc and Tc2 = Tscf is 
also present on the computed phase diagrams, e.g. on 
Fig.[T]3 for g> 0.65 VeV. 

The phase diagrams of LaGd and ThGd alloys, for 
which the free energy is given by Eq. (j76p with s = 7/2, 
are depicted in Figs. [5] and El The values of parame- 
ters exploited in these computations are collected in Ta- 
ble m Figures [5] and [3] show that the perturbative effect 
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FIG. 1: (Color online) Phase diagrams of (Lai_a;Cea;)Al2 for the values of M, g, S and GqPf given in Table|T]and under varying 
x: (a) X = 0.0010, (b) x — 0.0019, (c) x = 0.0028 and (d) x — 0.0057. The following convention was used in order to distinguish 
the regions corresponding to different phases: SC (red), D (blue), F (black) and the gapless phase NG (grey), where the 
smallest excitation energies from ground state disappear. 



Alloy 


X 


c[%] 


M 


g[VeV] 


S [eV] 


GapF 


(LaCe)Al2 


0.0010 


- 


1 


0.10 


0.01 


0.2610 


(LaCe)Al2 


0.0019 


- 


4 


0.189 


0.01 


0.2515 


(LaCe)Al2 


0.0028 


- 


7 


0.19 


0.01 


0.2435 


(LaCe)Al2 


0.0057 


- 


8 


0.23 


0.01 


0.2250 


LaGd 


— 


0.20 


5 


0.16 


0.01 


0.2825 


LaGd 


- 


0.30 


10 


0.18 


0.01 


0.2750 


LaGd 


- 


0.40 


15 


0.20 


0.01 


0.2670 


ThGd 


- 


0.10 


5 


0.08 


0.01 


0.2090 


ThGd 


- 


0.20 


8 


0.09 


0.01 


0.1910 



TABLE L Parameter values. 



of Gd is larger than that of Ce impurities, viz., the area 
of D phase is smaller than in Fig. [1] Furthermore, su- 
perconductivity is already suppressed at smaller values 
of g. This complies with experimental observations by 
Matthias et al., who showed that indeed the depression 
of superconductivity of doped lanthanum increases with 
the spin of the rare earth ionsi. 

The TscF temperature increases almost linearly with 
the value of g. One also observes decrease of Tpo with 
g and increase of Tpp with g. Furthermore, rapid dis- 
appearance of D phase is observed for critical values of 
9- 



The intermediate phase D in which superconductivity 
coexists with ferromagnetism of impurities is present on 
each phase diagram depicted in Figs. [TJ [5] and |31 The 
first experimental suggestions concerning the presence of 
coexistence phase in superconducting alloys was made by 
Matthias et ali. This hypothesis was later confirmed in 
Cei_^Gda;Ru2^ and Yi_2;Gda;0s2^. However, the basic 
question, as to whether superconductivity and magnetic 
order occurs in the same volume element was still open. 
Wilhelm and Hillenbrand measured phase diagrams of 
Cei_a;Tb2:Ru2, which contain the coexistence phase'^'^=. 
As temperature was lowered, one observed the following 
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FIG. 2: (Color online) Phase diagrams of LaGd for the values 
of M, (J, & and GopF given in Table U and under varying Gd 
concentration : (a) c = 0.10 % Gd, (b) c = 0.20 % Gd and 
(c) c = 0.30 % Gd. 



phase transitions, depending on the value of impurity 
concentration c: 

1. A phase transition P -^ SC, which occurs for suf- 
ficiently small c. 

2. For intermediate values of c the system undergoes 
a phase transition P — > SC and then from SC to 
D phase. 

3. Two phase transitions occur for sufficiently large 
values of c: P -^ F and then F ^ D. 

4. For some special value of c the system undergoes a 
phase transition P ^ D. 




T[K] 



b) 




T[K] 

FIG. 3: (Color online) Phase diagrams of ThGd for the values 
of M, g, 6 and GopF given in Table |l] and under varying Gd 
concentration : (a) c = 0.10 % Gd, (b) c = 0.20 % Gd. 



The phase diagrams depicted in Figs.[Tl[2]and[3]reveal the 
transitions 1-3. In particular, one phase transition P — >■ 
SC (variant no. 1) is visible on every diagram depicted in 
Fig. [1] for g < 0.05 VeV, which states the weak coupling 
of Cooper pairs to impurity ions. The options 2 and 3 are 
clearly seen on every diagram. However, variant no. 3 is 
most distinct in Fig. [1^ for g « 0.9 veV. The computed 
phase diagrams do not exhibit the fourth transition only, 
although, Fig.[T^ shows that the phase transition P — ^ D 
may appear for sufficiently large values of g, CqPf and 
small M and c, since in this case the D region is relatively 
large. 

The specific heat measurements of Cei_a;Tb2;Ru2 
revealed a short-range magnetic order in this com- 
pound, which is typical of spin-glasses^. The coex- 
istence of superconductivity with long-range ferromag- 
netic order was discovered a few years later, e.g. in 



(Eri_^H0:r)Rh4B, 



,52-54 



Phase diagrams of these com- 



pounds show the following phase transitions: 

1. P^ SC for X sa 0.1 

2. P ^ SC ^ D ioT X ^ 0.3. 

3. P ^ SC ^ F ^ D ioT X R:: 0.4. 
A. P ^ F ^ D ioT X ^ 0.9. 

5. P -i' F ioT X Ki 0.95. 
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Theoretical phase diagrams depicted in Fig. [THS] do not 
reproduce the third scenario only. 

The intermediate phase, where superconductivity co- 
exists with long-range antiferromagnetic order was dis- 
covered in RMoeSes (R = Gd, Tb and Er)^^, RRh4B4 
(R = Nd, Sm and Tm>33 and RMogSg (R = Gd, Tb, 
Dy and Er)-^. Upper critical field Hc2 in most of these 
compounds decreases below Neel temperature Tn. Ac- 
cordingly, the superconducting state is perturbed by 
antiferromagnetically aligned impurity spins. However, 
a rapid increase of Hc2 below Tn was also observed, e.g. 
in SmRh4B4, GdMoeSg, ThMogSg, Sni-^Eu^rMogSg, 
Pbi-a^EuxMoeSs, hai.2-x^nxMoeSi^^^, implying en- 
hancement of superconductivity by antiferromagnetism 
of impurities. 

Fischer et ali^ have pointed out that this boost of 
superconductivity may be due to the Jaccarino-Peter 
effecli^, in which superconductivity is induced by apply- 
ing a high magnetic field %. This effect can arise in a 
type-II superconductor, in which the impurity magnetic 
moments are antiferromagnetically coupled to conduc- 
tion electrons. This interaction generates an exchange 
field "H J , which acts on the spins of conduction electrons 
equivalently to an applied magnetic field, viz., breaks the 
Copper pairs. However, the negative sign of the cou- 
pling between the magnetic moments and the conduc- 
tion fermions spins, determines the direction of "Hj to be 
opposite to that of %. Thus, an applied magnetic field 
will be compensated by an exchange field, since the net 
magnetic field Ht is given \>y % — \'H}\. A given com- 
pound displays superconducting properties as long as the 
following relation holds: 



-T~Lp < ^T ^ ^p, 



where 



•H. 



PF 



XP - Xsc 



(77) 



(78) 



Xp and xsc denote the susceptibility of the normal 
and superconducting state respectively. "Hp, defined by 
Eq. ([75)) is the Chandrasekhar-Clogston limiting param- 
agnetic field^i^. 

The Jaccarino-Peter mechanism was observed exper- 
imentally in Euo.75Sno.25M06S7.2Seo. 3^. In the low- 
temperature scale and with increasing applied magnetic 
field four phase transitions were recognized: SC — > P — >■ 
SC — > P. This effect can explain the reappearance of su- 
perconductivity for sufficiently large values of g on theo- 
retical phase diagrams (Figs. Ill H and O. For T « O.OIK 
it is firmly seen that the system undergoes the following 
transitions: SC -^ D ^ F ^ SC, which can be inter- 
preted as SC — > P — > SC, since the given compound 
displays features characteristic of superconductors and 
magnetically ordered systems in a Z? phase and in the F 
phase superconductivity is suppressed. 

In conclusion it is worth pointing out that the interplay 
between superconductivity and magnetism is believed to 



be a possible mechanism of high-Tc superconductivity^ -, 
since the undoped state of cuprate superconductors is 
a strongly insulating antiferromagnet. The existence of 
such a parent correlated insulator is viewed to be an es- 
sential feature of high temperature superconductivity. 



VIII. CRITICAL TEMPERATURE 

According to section|Vl Eq. ([5^ for the solution {A 7^ 
0,1^ — 0} reduces to the BCS gap equation 



1 r Abcs 
Abcs = t^GqPf / -^ 

^ J-S £/BCS 



tanh(i/3PBCs)clC, (79) 



E, 



BCS 



- ./<2 



e 



^BCS- 



-.(BCS) 



The transition temperature Tc in BCS theory, is de- 

fined as the boundary of the region beyond which there 

is no real, positive Abcs satisfying ([79]) . Below Tc 
the solution Abcs 7^ minimizes the free energy and the 
system is in superconducting phase. Therefore, Tc 
can be obtained from Eq. ([75| with Abcs = Oj which 
yields^: 



T^^'^^^^l.USe^pl-iGoPv)-']. 



(80) 



It should be possible to estimate the change in Tc , 
since the density of states enters exponentially in (|80p . 
However, significant deviations from Eq. (IMl) were ob- 
served experimentally for a number of superconductors 
containing magnetic impurities. This inadequacy of 
Eq. ([50)1 is most distinct for large values of impurity con- 
centration. BCS theory is therefore incapable to describe 
the superconducting alloys. 

The results obtained in Sec. IVllI show that phase tran- 
sition to superconducting state in a superconductor with 
magnetic impurities, depending on the value of magnetic 
coupling constant, can be 1st or 2nd order. Next two 
subsections are concerned with computation of the tran- 
sition temperature Tc on the grounds of Eqs. (|69p . (170)) 
and (|76[) . These calculations strongly depend on the or- 
der of the transition. 



A. Second order phase transition 

Expression for transition temperature Tc for 2nd order 
phase transition can be computed analogously as in BCS 
theory. It suffices to put A = in Eqs. (HH) and (|70)). 
Thus, one obtains the following set of equations for Tc = 

l/A:/3c, 



GqPf 



d^ 



sinh(/3c|CI) 



s ICI cosh(/3c|ei) + co8h{gpcMf^''\,yc)) ' 



(81) 
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eg shMgp,Mf!f\v,)) 

M cosh{gp,Mf!f\v,)) + cosh(/3e|ei) 



+ /i^^(^c), 



(82) 

where Vc — i^(/3c), s = 1/2, 7/2. Numerical analysis show 
that in the low-temperature scale z^c is almost indepen- 
dent in T, viz., 



for spin 1/2 impurities: 



0, for X < 0.0010, 



1/(0) = ^, for a; > 0.0019 



for spin 7/2 impurities: 



^c«.(0) = f. 



This result complies with experimental data, showing 
that the perturbative effect of impurities on supercon- 
ductivity is an increasing function of their spin. Fur- 
thermore, the reduction of Tc by adding a small amount 
of spin 1/2 impurities can be satisfactorily described by 
Eq. (IMl) by adjusting the value of Gq/Of, since for i^c = 
Eqs. dH]) and dH]) reduce to T^^'^^\ Accordingly, the set 
of Eqs. (|8ip. (|5^ are solved under the assumption that 

w — £9 

i^c — M- 

The solution for Tc resulting from Eqs. ([81]), ^ for 

s = 1/2, 7/2 and under varying impurity concentration 
c is depicted in Figs. SI O In general, the values of the 
parameters M, g, S and GqPf, depend on the impurity 
concentration c. However, to minimize the number of 
adjustable parameters, the values of Af , g, 6 and GoPf 
have been kept constant. 

Figures |4j [5] show how variation of c affects Tc (c) . For 
certain values of c (e.g. c G (0.0225, 0.0275) for curve 
no. 1 on Fig. SJ Eq. (JST]) with Vc ~ ^^(0) possesses two so- 
lutions: Tci, rc2. Thus, superconductivity is suppressed 



below Tc2- The dependence of Tc on c shown in Figs. lU [5] 
does not indicate the existence of a third transition tem- 
perature Tc3. This, can be a consequence of the applied 
approximation in which M, g, S i GqPf do not depend on 
c. However, the dependence of Tc on impurity concen- 
tration c shown in Fig. S] is typical of the Lai_2;Cea;Al2 
alloyl^. 

The destructive influence of impurities on supercon- 
ductivity is analogous to the effect of an external mag- 
netic field H applied to a superconducting compound. 
Sarma§2 obtained a numerical solution for Tc{H) of the 
system described by a Hamiltonian Hs = -ffBCS+Z^B'Hcrz. 
His result for Tc('H) agrees qualitatively with Tc(c) 
graphs depicted in Figs. S] and [51 since the expression 
for Tciji) obtained in Ref. |63| is of the similar form to 
Eq. dHI]) with gMf^'^Vc) replacing ^bH. 



B. The first order phase transition 

In the case of 1st order phase transition, the assump- 
tion that the gap parameter A vanishes at the transition 
temperature does not hold. According to the results ob- 
tained in the previous section, one can expect that Tc 
possesses three solutions (Td < Tc2 < Tcs) for certain 
values of g and c. These solutions can be determined 
numerically from the following equations 



T 



cl 



Fp — Fsc — 0, 



Tr, 



Fsc-F^^O, ^^D,F, 



T 



c3 



Ff -Fn = 0. 



(83) 



(84) 



(85) 



The existence of Tc3 depends on the type of phase tran- 
sition occuring at Tc2. If the system undergoes a phase 




0,000 0,005 0,010 0,015 0,020 0,025 0,030 




0,000 0,002 0,004 0,006 0, 



0,010 0,012 0,014 



FIG. 4: (Color online) The superconducting transition tem- 
perature Tc under varying impurity concentration c for su- 
perconducting alloy, containing 1/2 impurities. The values of 
the parameters g, GopF and Af are given in table HIl 



FIG. 5: (Color online) The superconducting transition tem- 
perature Tc under varying impurity concentration c for su- 
perconducting alloy, containing 7/2 impurities. The values of 
the parameters g, GopF and M are given in table HIl 
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Impurity spin 


Curve no. 


M 


5[VW] 


5 [eV] 


GopF 


1/2 


1 
2 

3 

4 


1 

4 

7 
8 


0.10 
0.189 
0.19 
0.23 


0.01 
0.01 
0.01 
0.01 


0.2610 
0.2515 
0.2435 
0.2250 


7/2 


1 
2 
3 


5 
10 
15 


0.16 
0.18 
0.20 


0.01 
0.01 
0.01 


0.2825 
0.2750 
0.2670 



TABLE II: Parameter values. 



m 



qn 



^ 



0,000 0,001 0,002 0,003 0,004 0,005 0,006 0,007 



0,5 

0.0 



^^ 



*r 



0,000 0,001 0,002 0,003 0.004 0,005 0,006 



FIG. 6: (Color online) The superconducting transition tem- 
perature vs Ce impurity concentration x for (Lai_a;Cea;)Al2. 
Squares denote theoretical values of Tc obtained from 
Eqs. (I83p - (|85|l . Rhombs are experimental points from Ref.ilSl. 
The star at the point x — 0.0067 denotes the experimentally 
estimated turning point in Tc, above which superconductivity 
is destroyed. 



transition to ferromagnetic phase at Tc2, then T^a > 
for certain values of g (Sec. IVIip . If Tc2 — Tscd, then 
Tc3 — and the system does not reenter the supercon- 
ducting phase {SC or D). 

In order to obtain a reliable comparison of the tran- 
sition temperature resulting from Eqs. ([55)) ~(|55 |) to ex- 
perimental data, we let the parameters Gopp, g and 
M vary with impurity concentration. The parameters 
are then adjusted independently for each experimen- 
tal point. The best fitting of Tdc) to experiment for 
Lai_2;Ce2,Al2, Lai_2;Gdj;Al2 and (Lao.8-a;Yo.2o)Cea; is 
plotted in Figs. El [7] and [H In all cases very good quan- 
titative agreement with experiment was found. 

Fig. [S] shows two phase transitions for 

X e (0.0064, 0.0067). The first phase transition 
[P -^ SC) occurs at Td and the second transition 
(back to normal state) appears at Tc2- The reentrant 
superconductivity (the transitions SC — >■ P — >■ SC) 



FIG. 7: (Color online) The superconducting transition tem- 
perature vs Gd impurity concentration x for (Lai_j;Gda;)Al2. 
Squares denote theoretical values of Tc obtained from 
Eqs. (I83l) - (|85p . Rhombs are experimental points from Ref . l6a. 
The value of critical concentration a^cr, where Tc = is 0.59 
at. % Gd. 



is clearly visible in (Lao 8-2; Yo 2o)Ce2; for x = 0.0085 
(Fig.©. 

The theoretical curve Tc{x) resulting from AG the- 
ory also yields very good quantitative agreement with 
experimental data for Lai_a;Gd2:Al2. However, sig- 
nificant deviations from this theory were observed in 
(Lai_a;Cea;)Al2 and (Lao.8-a;Yo.2o)Cea;, since AG pre- 
dicts only single- valued solutions for the superconducting 
transition temperature^^. The model studied here gives 
considerably better results for Tc than previous theories, 
although analytic expression for Tc cannot be obtained 
for 1st order phase transitions. 



IX. CONCLUDING REMARKS 

We have shown that the thermodynamics of a BCS 
superconductor, perturbed by a reduced s-d interaction 
is solvable if the impurity density is sufficiently small. 
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FIG. 8: (Color online) The superconducting transition temperature vs Ce impurity concentration x for (Lao.8-a;Yo.2o)Cea:. 
Squares denote theoretical values of Tc obtained from Eqs. (|83p - (|85|l . Rhombs are experimental points from Ref. Ha . 
For X — 0.0085 there are three Tc solutions, which demonstrate reentrant superconductivity observed experimentally in 
(Lao.8-a;Yo.2o)Cea; at this value of a; (inset). 



The essential properties of dilute magnetic supercon- 
ducting alloys have been demonstrated: the decrease 
of superconducting transition temperature Tc with in- 
creasing impurity concentration, the existence of reen- 
trant and gapless superconductivity and the presence 
of Jaccarino-Peter compensation related to the mag- 
netic field induced superconductivity. Good quantitative 
agreement of the resulting dependence of T^ on impurity 
concentration was demonstrated for several supercon- 
ducting compounds: Lai_a;Cea;Al2, Lai_2;Gd2;Al2 and 
(La.Q,s-x^a.2o)Cex- The theory presented provides better 



agreement with experiment than earlier AG, MHZ the- 
ories and their various extensions and refinements given 
in Refs. [l^-ilH. 

We have also performed a fit of specific-heat and crit- 
ical field curves for Lai_a;Cea;Al2, LaCe, LaGd, and 
ThGd. We shall report on this study elsewhere^. 

These investigations will be extended to include the ef- 
fect of a general s-d exchange interaction T4_(j and a BCS- 
type attraction between Cooper pairs ¥41^^— in order to 
study other properties of superconducting alloys. 
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